Neutral skyrmion configurations in the low-energy effective theory of spinor 

condensate ferromagnets 
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We study the low-energy effective theory of spinor condensate ferromagnets for the superfiuid 
velocity and magnetization degrees of freedom. This effective theory describes the competition 
between spin stiffness and a long-ranged interaction between skyrmions, topological objects familiar 
from the theory of ordinary ferromagnets. We find exact solutions to the non-linear equations 
of motion describing neutral configurations of skyrmions and anti-skyrmions. These analytical 
solutions provide a simple physical picture for the origin of crystalline magnetic order in spinor 
condensate ferromagnets with dipolar interactions. We also point out the connections to effective 
theories for quantum Hall ferromagnets. 



I. INTRODUCTION 

For systems with broken symmetries, low-energy ef- 
fective theories provide a simple and powerful tool for 
describing the relevant physics. The focus is on the long- 
wavelength Goldstone modes which emerge from the mi- 
croscopic degrees of freedom. This often reveals the con- 
nections between seemingly unrelated systems that share 
the same pattern of symmetry breaking. The most well- 
known example is that of the complex scalar field used 
in the low-energy theories of bosonic superfluids [3, 13, 
fermionic superconductors Q, and XY spin systems j4|. 

Effective theories also bring to light subtle topological 
effects that may become important at low energies. Gold- 
stone modes often carry a non-trivial topology which can 
give rise to the appearance of topological defects. For ex- 
ample, in two spatial dimensions Kosterlitz and Thouless 
pointed out the crucial role of vortices for the complex 
scalar field 0, [|| . Vortices are point- like topological de- 
fects around which the phase of the complex scalar field 
winds by an integer multiple of 2ir. 

Ultracold atomic systems provide an ideal testing 
ground for low-energy effective theories. The microscopic 
degrees of freedom are well-isolated from the environment 
experimentally and well-understood theoretically. The 
challenge is in describing how these microscopic degrees 
of freedom organize at low energies in the presence of 
non-trivial interactions. When only one internal hyper- 
fine level is important, the phenomenon of scalar Bose- 
Einstein condensation is given by the low-energy theory 
of the complex scalar field [6|-|8[ • A series of ground break- 
ing experiments have observed vortex lattices in rotating 
condensates 0, [l(| as well as evidence for the role of vor- 
tices in the equilibrium Kosterlitz-Thouless transition for 
two-dimensional condensates 

For ultracold atoms with a complex internal level struc- 
ture, there are various patterns of symmetry breaking. 
This lead to rich possibilities and challenges for low- 
energy effective descriptions. Spinor condensate ferro- 
magnets are one such system which has seen a number of 
important experimental advancements (see Ref. ,12]). 
This includes the development of optical dipole traps 



used for preparation [13[ and phase-contrast imaging 
used for detection 14| in S = 1 87 Rb. In addition to the 
phase degree of freedom familiar from single-component 
condensates, the magnetization naturally arises as a de- 
scription of the low-energy spin degrees of freedom. A 
vector quantity sensitive to both the population and co- 
herences between the three hyperfine levels, the magne- 
tization can be directly imaged in experiments [ldj . 

One of the most striking observations in spinor con- 
densate ferromagnets is the spontaneous formation of 
crystalline magnetic order [la , Il6j . From an initial 
quasi-two-dimensional condensate prepared with a uni- 
form magnetization, a crystalline lattice of spin domains 
emerges spontaneously at sufficiently long times. The 
presence of a condensate with magnetization sponta- 
neously breaks global gauge invariance and spin rota- 
tional invariance. Additionally, crystalline order for the 
magnetization breaks real space translational and rota- 
tional symmetry. Previous works have pointed out the 
crucial role of dipolar interactions in driving dynami- 
cal instabilities within the uniform condensate towards 
states with crystalline order [l7i - fT9j [ . Numerical analysis 
of the full multi-component Gross-Pitaevskii equations 
suggest dipolar interactions can give rise to states with 
crystalline order [18|, l20( • 

In this and a companion paper (2l| . we take a com- 
plementary approach and focus on the low-energy effec- 
tive theory of two-dimensional spinor condensate ferro- 
magnets. This effective theory describes the interaction 
between the superfiuid velocity and the magnetization 
degrees of freedom. Previous work has derived the equa- 
tions of motion for this effective theory [2^, [H|. We 
extend this result to demonstrate how the Lagrangian 
for the effective theory can be written as a non-linear 
sigma model in terms of the magnetization alone. The 
effect of the superfiuid velocity is to induce a long-ranged 
interaction term between skyrmions, topological defects 
familiar from the theory of ordinary ferromagnets [24j . 
In contrast to the point-like topological defects of vor- 
tices, skyrmions describe extended magnetization tex- 
tures which carry a quantized topological charge. 

In the companion paper [2lj , we study how symmetry 
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groups containing combined real space translational, real 
space rotational, and spin space rotational symmetry op- 
erations can be used to classify possible crystalline mag- 
netic orders. Within each symmetry class, we find min- 
imal energy configurations describing non-trivial crys- 
talline configurations. 

The main purpose of this paper is to give a simple 
physical physical picture for the origin of crystalline or- 
der in terms of neutral configurations of skyrmions and 
anti-skyrmions. Long-ranged skyrmion interactions force 
magnetization configurations to have net neutral collec- 
tions of topological defects. We show this explicitly 
by finding exact analytical solutions for the non-linear 
equations of motion describing both localized collections 
of skyrmions and anti-skyrmions as well as extended 
skyrmion and anti-skyrmion stripes. 

Skyrmions have non-trivial spin configurations that 
spontaneously break translational and rotational invari- 
ance in real space. Proposed originally in high en ergy 
physics as a model for mesons and baryons (2|| l26j |. 
skyrmions have found applications in a number of diverse 
fields including quantum hall ferrogmagnets [27l [28| , and 
magnetically ordered crystals |29l432j . 

A neutral collection of such topological objects is 
able to take advantage of the dipolar interaction en- 
ergy without a large penalty in the skyrmion interac- 
tion energy. Since a neutral skyrmion configuration has 
the same topological number as the uniform magnet, its 
stability is not ensured by topology alone. However, 
the scale invariance of the skyrmion interaction energy 
rules out the most straightforward instability of bring- 
ing skyrmions and anti-skyrmions closer together. Essen- 
tially, the skyrmion interaction energy forces the charge 
densities of a skyrmion and anti-skyrmion pair to shrink 
as the distance between them shrinks so that the overall 
energy remains the same. 

The analytical solutions we find without dipolar inter- 
actions closely resemble the minimal energy configura- 
tions found numerically in the presence of dipolar inter- 
actions. The role of dipolar interactions can then be seen 
as stabilizing these non-trivial solutions of the effective 
theory. We point out that magnetic dipolar interactions 
are small. Thus it is a good starting point to find states 
which are static nontrivial solutions of the system with- 
out dipolar interactions. 

The effective theory of spinor condensate ferromag- 
nets is essentially identical to that of the quantum Hall 
ferromagnets [23, [Hj]. Although the microscopic de- 
grees of freedom are fermionic electrons, the magneti- 
zation order parameter is bosonic. The Coloumb inter- 
action between electrons then gives a contribution to the 
skyrmion interaction for the effective theory. The result- 
ing skyrmion interaction is qualitatively the same as the 
one for spinor condensate ferromagnets. This suggests 
the study of spinor condensate ferromagnets may have in- 
teresting connections to quantum hall ferromagnets and 
vice versa. 

The plan of this paper is as follows. In Sec. |H] we 



review the theory of ordinary ferromagnets and how 
skyrmion solutions arise from the non-linear Landau- 
Lifshitz equations of motion. We then proceed to re- 
view how these skyrmion solutions are used in the low- 
energy effective theory of quantum Hall ferromagnets in 
Sec. Mil Next we derive the low-energy effective the- 
ory for spinor condensate ferromagnets in Sec. IIVI In 
particular, we demonstrate how a long-ranged skyrmion 
interaction term (which also appears for quantum hall 
ferromagnets) arises from coupling of the magnetization 
to the superfluid velocity. 

In Sec. [Vj we discuss how to interpret the mathemati- 
cal structure of skyrmion solutions for ordinary ferromag- 
nets in terms of a separation of variables. This approach 
allows us to use find new exact solutions for the oridinary 
ferromagnet. More importantly, it also allows us to gen- 
eralize the skyrmion solutions of the ordinary ferromag- 
net to find analytical solutions for the spinor condensate 
ferromagnet. These latter solutions describe both neutral 
collections of localized skyrmions and anti-skyrmions as 
well as extended stripe configurations. Finally, we dis- 
cuss how the analytical solutions we find offer insight 
into quantum hall ferromagnets and spinor condensate 
ferromagnets with dipolar interactions in Sec. I VII 



II. SKYRMIONS IN FERROMAGNETS 

We begin by reviewing the theory of ordinary two- 
dimensional ferromagnets described by the following La- 
grangian, Hamiltonian, and Landau-Lifshitz equations of 
motion [34j 

C = -S J dtd 2 xA(h) -d t n- J dtH 

H = ^J d 2 xV{hf 

d t h = -h x V 2 n (1) 

where n is a three component real unit vector and A(h) 
is the unit monopole vector potential. The order param- 
eter h describes the magnetization and is a unit vector 
living on the sphere. Calculating the variation of the La- 
grangian to derive the Landau-Lifshitz equations [35l . I3H ] 
can be done by using Sn — Sw x h. This is consistent 
with the constraint |ri| = 1 since Sh ■ h = by construc- 
tion. The variation of the term J dtA(h) • dth is clearest 
when using its geometric interpretation as the area on 
the sphere swept by h. While we include time depen- 
dent terms for completeness, in this paper we will only 
consider static solutions 

In addition to the trivial uniform solution, there are 
non-trivial soliton solutions to the Landau-Lifshitz equa- 
tions called skyrmions. By parameterizing 

ft = [sin(a) cos(/3) sin(a) sin(/3) cos(a)] T (2) 




FIG. 1: A single localized skyrmion carrying 4-7T net skyrmion 
charge in the ordinary ferromagnet (left). Neutral configura- 
tion consisting of a localized skyrmion carrying +4-7T skyrmion 
charge in a negative background carrying — 4tv skyrmion 
charge in the spinor condensate ferromagnet (right). Notice 
the +z (— z) meron carrying +2tt (+2tt) net skyrmion charge 
at the origin (at infinity) for the ordinary ferromagnet are 
mapped to a skyrmion (anti-skyrmion) carrying +4tt (-47r) 
net skyrmion charge in the spinor condensate ferromagnet. 
Red (blue) background indicates positive (negative) skyrmion 
density q, black 2D arrows the superfluid velocity v for the 
spinor condensate ferromagnet, and shaded 3D arrows the 
magnetization h. 

we see a controls fi z , the z component of the magnetiza- 
tion while ft controls the canonically conjugate variable 
giving the orientation of n x , h y , the x, y components of 
the magnetization. 

The minimal energy solutions within each topological 
sector are called skyrmions and can be written in the 
form [H 

tan(a/2)e^ = exp[/(x + iy)] = exp[it(x, y) + iv(x, y)] 

(3) 

where f(z) is a holomorphic function of z with real part 
u(x,y) and imaginary part v(x,y). The function f(z) ~ 
nlog(z — zo) can have logarithmic singularities [52|j. Since 
j3 has to be 2ir periodic, the residues of these singularities 
must be integers. This implies exp[/(x + iy)] can only 
have zeros or poles. The resulting spin configuration has 
n z = +1 (n z = —1) at zeros (poles) of exp[/(z)] while h x , 
fiy wind anti-clockwise (clockwise) along a path circling 
the zeros (poles) in a anti-clockwise direction. 

These solutions describe topological defects of ordinary 
ferromagnets. For fixed boundary conditions, smooth, 
finite energy configurations are separated into distinct 
classes characterized by a quantized topological invariant 
J d xq(x) = AttN where N is the number of times n 
covers the sphere. Here 

q = £ F n • V M n x V„n (4) 

is the skyrmion density which is positive for f(z) holo- 
morphic. From here on, lower Greek indices (upper Ro- 
man indices) refer to real space (order parameter) com- 
ponents. 

Consider the single skyrmion solution shown in the left 
of Fig. [T] It corresponds to f(z) = log(z), carries net 



FIG. 2: Unit cell for a lattice of localized skyrmions carrying 
87r net skyrmion charge per unit cell in the ordinary ferromag- 
net (left). Unit cell for a neutral configuration carrying zero 
net skyrmion charge per unit cell in the spinor condensate 
ferromagnet (right). Notice +z ( — z) merons are mapped to 
skyrmions (anti-skyrmions). Red (blue) background indicates 
positive (negative) skyrmion density q, black 2D arrows the 
superfluid velocity v, and shaded 3D arrows the magnetiza- 
tion n. 



skyrmion charge 47r, and can be decomposed into a +z 
meron at the origin and a — z meron at infinity. Essen- 
tially, a meron can be thought of as half of a skyrmion 
and characterized by two signed quantities: the direction 
of the magnetization at the core and the orientation of 
the winding away from the core. The sign of the skyrmion 
density is the product of the sign of these two quantities. 
For example, the +z meron at the origin to the left of 
Fig. [1] carries net skymrion charge 2ir. 

For spinor condensate ferromagnets, scale invariance 
of the skyrmion interaction term guarantees stability 
against trivial rescaling. However, this may change if 
we introudce a short distance cutoff or quantum fluctua- 
tions. 

Notice the Lagrangian in Eq. [1] has translational, ro- 
tational, and scale invariance in real space as well as ro- 
tational invariance in spin space. The single skyrmion 
solution spontaneously breaks all of these symmetries. 
However, the action is invariant for f(z) = \og[fo(z — zq)] 
with complex constants zq and Jo . Taking solutions with 
20 7^ corresponds to spatially translating the solution 
with zo = 0, solutions with Re[/o] 7^ 1 correponds to 
spatially rescaling the the solution with /o = 1, and so- 
lutions with Im[/o] 7^ corresponds to rotation of the 
/o = 1 solution. In addition, the action is the same for 
On where O is a rotation matrix describing solutions re- 
lated by spin space transformations. 

Collections of multiple skyrmions arc described by 
f(z) having multiple singularities. A square lattice of 
skyrmions with net skyrmion charge 8ir per unit cell is 
shown in the left of Fig. O It can also be decomposed 
into a collection of two +z merons and two — z, merons 
arranged antiferromagnetically. In addition, there are 
solutions tan(a/2)e 1 ^ = exp[/(ai — iy)] with f(z) anti- 
holomorphic characterized by x, y components winding 
in the opposite direction compared to holomorphic solu- 
tions and q negative. 
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III. QUANTUM HALL FERROMAGNETS 



IV. SPINOR CONDENSATE FERROMAGNETS 



After describing the structure of skyrmion solutions 
for ordinary ferromagnets, we now briefly review how 
they are used in the study of quantum Hall ferromagnets 
[27l l33j . For the low-energy effective theory of quantum 
Hall ferromagnets, spin 1/2 electrons in a magnetic field 
are described by a two component Chern-Simons com- 
posite boson theory. The bosons couple to the phys- 
ical gauge field due to the magnetic field and a ficti- 
tious Chcrn-Simons gauge field which attaches appro- 
priate flux quanta to change from bosonic to fermionic 
statistics. At appropriate filling fractions, the physical 
and Chern-Simons fluxes cancel and the resulting quan- 
tum Hall plateau is described as a condensate of compos- 
ite bosons. 

The Lagrangian and Hamiltonian for the resulting ef- 
fective theory is given by 



C = 

n = 



dtd 2 xA(n) ■ dtfi 



dm 



d 2 3 



!(v«) 2 



-B ■ n 



d 2 xd 2 y[q(x)~q]G{x-y)[q{y)~q] (5) 



where B is the magnetic field giving rise to a linear Zee- 
man shift, q(x) is the skyrmion density and q = J d 2 xq(x) 
is the net skyrmion charge. 

The coupling of the bosons to the Chern-Simons field 
gives rise to the skyrmion interaction term. At long dis- 
tances, the Chern-Simons field ties the skyrmion density 
to the deviation of the physical electron density from its 
background value. In fact skyrmions are lowest energy 
quasiparticles at a filling factor of one [13, H3]- Thus, the 
skyrmion density inherits the Coloumb interaction for the 
electron density. Away from quantum Hall plateaus, the 
background value q is non-zero and the singular \x— 
behavior in G{x — y) forces finite energy configurations 
to have non-zero net skyrmion charge. 

Thus the focus is on static configurations carrying net 
skyrmion charge. The skyrmion solutions of the ordinary 
ferromagnet have non-zero q and provide a good qualita- 
tive description of quantum Hall ferromagnets away from 
quantum Hall plateaus. For example, a sqaure lattice of 
skyrmions as shown in Fig. [5] carries a net charge of 
8ir per unit cell and can be used as a starting point for 
studying configurations carrying finite q. Although such 
solutions take into account the spin stiffness and the long- 
ranged divergence of the skyrmion interaction, quantita- 
tive results require additional analysis arising from the 
detailed form of skyrmion interaction and the linear Zee- 
man shift. This usually involves numerical minimization 
of the Hamiltonian in Eq. [3] 



Having reviewed known results on ordinary and quan- 
tum Hall ferromagnets, we now consider spin S spinor 
condensate ferromagnets described by the microscopic 
condensate wavef unction a 2S + 1 complex vector 
[4ll [42[. The microscopic Gross-Pitaevskii Lagrangian 
is given by 



C = I dti&dt* 
H = I d 2 x 



dm 



dtHt 



1 

2m 1 



■g.(WF¥) 



(6) 



where F are spin matrices and go > gives the spin- 
independent contact interaction strength while g s < 
gives the spin-dependent contact interaction strength 
which favors finite magnetization. Here, Us denotes 
additional spin dependent interactions such as the 
quadratic Zeeman shift and dipolar interactions. 

As was the case for quantum Hall ferromagnets, we 
expect a simpler description to emerge at low energies. 
The resulting low-energy effective theory should only in- 
volve the condensate phase (f> and local magnetization 
n which describe the order parameters of the system. 
Previous work has shown this can be done at the level 
of the equations of motion 22] . In this section, we ex- 
tend this result to derive the Lagrangian and Hamilto- 
nian for the effective theory solely in terms of the magne- 
tization. However, the skyrmion density acts as a source 
of voriticity for the superfluid velocity. Thus, the ef- 
fect of the superfluid phase is to induce a logarithmic 
vortex-vortex interaction between skyrmions. The re- 
sulting non-linear sigma model is essentially identical to 
that of the quantum Hall ferromagnet but with skyrmion 
interaction G(x — y) ~ log(x — y) having logarithmic be- 
havior instead of \x — behavior. 

We begin by considering energies below the scale of 
spin-independent go and ferromagnetic spin-dependent g s 
contact interactions. The condensate has fixed density 
= p and fully polarized magnetization ^ = 
Spn where F are spin matrices. The states that satisfy 
these constraints are parameterized solely in terms of the 
low-energy degrees of freedom 



Ftpfi = Spi>fi 



(7) 



where </> descrribes the phase of the condensate and tpn is 
a fully polarized unit spinor with h describing the orien- 
tation of the magnetization. Although (f> is not directly 
observable, the superfluid velocity is a physical quantity 



(8) 



which has contributions from both <f> and ipn- 

In this paper, we are primarily interested in the com- 
petition between the spin stiffness and superfluid kinetic 



5 



energy. In the companion paper [21], we address the ef- 
fect of magnetic dipolar interactions. From here on, we 
consider the case Hs — 0. From the Gross-Pitaevskii 
Lagrangian in Eq. [51 the Berry's phase term becomes 
jxj/t^vEr = —pQ t( p— SpA{h) ■ dth while the kinetic energy 
term is given by |V\&| 2 = Sp/2(Vn) 2 +pv 2 and the inter- 
action terms give constants. This gives the Lagrangian 
and Hamiltonian as 



C = -S dtd z xA(h) ■ d t h 



dm 



H = d 2 x 



f(Vn) 2 + iv 2 



(9) 



where we take p = m = 1 for simplicity. Notice for fixed 
p, the d t (/) term is a total derivative which we exclude. 
Compared to the Lagrangian describing ordinary ferro- 
magnets in Eq. [T] there is an additional superfluid kinetic 
energy term v^v^. 

The global phase cf> enters the Lagrangian quadrati- 
cally and only through v. The equation of motion for 
cj> gives VpV^ = implying the superfluid velocity is di- 
vergenceless. This follows from v describing transport of 
the density p, a conserved quantity which is locally fixed 
at low energies due to the spin-independent contact in- 
teraction. This implies that v^v^ only depends on the 
divergenceless part of v. In momentum space, this is 
v Al (+fc)[(5 py — k^k v /k 2 \v u (—k) which can also be written 
as F f , l ,(+k)F llv (-k)/2k 2 . 

Here we have introduce the analog of the field strength 
tensor F M „ = V^v^-V^v^, a local quantity that depends 
only on the divergenceless part of v. In two dimensions, 
there is only one non-zero component to F^ u . From Eq. 
[5] this is given by the skyrmion density F xy = —F yx = 
Sq. Here we assume that the condensate phase <f> does 
not contribute to F xy through vortex-like singularities. 
This is valid because the vortex core energy is large. 

The above results give v solely in terms of h as 



0, 



Cuv'VuVv = Sq 



(10) 



with q the skyrmion density. Gradients in the order pa- 
rameter n arise in part from phase gradients in the con- 
densate wavefunction \&. Topologically non-trivial mag- 
netization configurations can thus give rise to vorticity 
described by a non-zero curl e^V^v,, ^ 0. 

By introducing the two-dimensional logarithmic 
Green's function — V 2 G(x) = S(x) we can write the 
superfluid kinetic energy F^(+k)F^(— k)/2k 2 in real 
space and obtain 

C = -S J dtd 2 xA(h) -d t h~ J dm 

H = | J d 2 x(Vh) 2 + ^-J d 2 xd 2 yq(x)G(x - y)q(y) 

(11) 



with the corresponding equations of motion given by 



1 



(d t + v^V^) h = - fix V z h 



(12) 



along with Eq. [10] for the superfluid velocity solved by 

v M = Se^V u <f>, -V 2 $ = q (13) 

where $>(x) = J d 2 xG(x — y)q(y) has the interpretation of 
the two-dimensional Coloumb potential associated with 

Compared to the Landau-Lifshitz equations describing 
ordinary ferromagnets in Eq. [IJ the replacement dt — > 
dt + v /iV M describes the advection of the magnetization 
by the superfluid velocity [22| . This advective term arise 
from variation of the superfluid kinetic energy term in Eq. 
[9] or equivalently from the skyrmion interaction term in 
Eq. [TT] Recall that we include the time dependence for 
completeness and focus only on static solutions. 

The skyrmion density q gives the vorticity for the 
superfluid velocity v. Thus, the second term in the 
Hamiltonian above gives the pairwise logarithmic inter- 
action energy between vortices. In the thermodynamic 
limit, the logarithmic divergence of G(x — y) at large dis- 
tances forces finite energy configurations to have zero net 
skyrmion density J d 2 xq(x) = 0. 

Notice Eq. [11] for spinor condensate and Eq. [5] for 
quantum hall ferromagnets have the same form. Al- 
though G(x — y) behaves as log \x — y\ for the former and 
| — 2/| 1 for the latter, both give singular contributions at 
long wavelengths. However, the important absence of a 
finite backgroud value q in the skyrmion interaction im- 
plies configurations for spinor condensate ferromagnets 
must have zero net skyrmion charge. 



V. EXACT SOLUTIONS WITH NEUTRAL 
SKYRMION CHARGE 



In quantum hall systems density deviations from the 
incompressible state cause skyrmions. Density fluctua- 
tions with zero net average (such as impurities) cause 
spin textures with zero net skyrmion number. States 
with non-zero skyrmion number occur away from quan- 
tum hall plateaus. Thus the analytical skyrmion solu- 
tions carrying net charge for the ordinary ferromagnet 
offered insight into more complicated case of quantum 
Hall ferromagnets away from the quantum Hall plateau. 
For spinor condensate ferromangets, we will show how 
net neutral solutions with skyrmions and anti-skyrmions 
without dipolar interactions offer insight into the more 
complicated case with dipolar interactions. 

In this section, we find exact analytical solutions 
for spinor condensate ferromagnets with logarithmic 
skyrmion interactions in the absence of dipolar inter- 
actions. We study the effect of including dipolar in- 
teractions numeri cally after a symmetry analysis in the 
companion paper [21(. The exact solutions we find here 
greatly resemble the numerical solutions in the compan- 
ion paper. As we discuss in Sec. I VII the interpretation 
of the exact solutions in terms of neutral collections of 
skyrmions and anti-skyrmions offers physical insight into 
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the more complicated numerical solutions of the compan- 
ion paper. 

To find exact solutions with zero net skyrmion charge, 
it is vital to include to effect of the skyrmion interac- 
tion term. Recall it is the long wavelength divergence 
of this term that forces configurations to have zero net 
skyrmion charge. Although this cannot be done exactly 
for \x — interactions as in quantum hall ferromag- 
nets, it is possible for logarithmic interactions as in spinor 
condensate ferromagnets. Physically, this is because the 
logarithmic interaction arises solely from the superfluid 
kinetic energy which is scale invariant just like the spin 
stiffness term. 

We begin with the parameterization of n in Eqs. [21 02 
used in the skyrmion solutions of the ordinary ferromag- 
net. Notice a, (3 provide a set of orthogonal coordinates 
for the sphere describing the order parameter space of 
h. For f(x + iy) — u{x,y) +iv(x,y) holomorphic, u(x,y) 
and v(x, y) provide a set of orthogonal coordinates for the 
plane describing real space. So for ordinary ferromag- 
nets, skyrmion solutions are given by a — 2tan _1 (e u ) 
and P = v (see also Eq. [3]), which can be understood as 
a separation of variables. 

Notice a(u) is a function of u only while P{v) is a 
function of v only. Each orthogonal coordinate of the 
order parameter space a, /? is a function of only one or- 
thogonal coordinate of real space u, v, respectively. The 
reason why using u and v as coordinates is tractable 
is because they satisfy the Cauchy-Riemann equations 
d x u = +d y v, d y u = —d x v. In particular, this implies 
Vu • Vw = meaning countour lines of constant u are 
perpendicular to countour lines of constant v as required 
for orthogonal coordinates. In addition, both V 2 u = 
and V 2 v = satisfy Laplace's equation. The above two 
identities simplify expressions involving V 2 which arise 
in the equations of motion. In particular, when changing 
variables from (x,y) to (u, v), the Laplacian retains its 
form <9 2 + <9 2 oc <9 2 + 9 2 . 

An alternative interpretation of the above separation 
of variables is as follows. Given an arbitrary configura- 
tion for h, consider the contour lines of constant h z , the z 
component of the magnetization. For contour lines with 
n z 7^ ±1 that form closed curves, consider the winding 
number of h x + ih y = sin(a)e 1 ^. For smooth configura- 
tions, this is a quantized integer that cannot change be- 
tween neighboring contours which do not cross n z =/= ±1. 
This implies the winding number is constant in regions 
between contours with h z ^ ±1. Label different contours 
of h z by u and the label coordinate along each contour 
by v. 

In order to have a non-zero winding number, /3 must 
have some dependence on v and the minimal one is p cx v. 
In principal, p can also depend on u and have some non- 
monotonic dependence on v, but linear dependence is the 
smoothest one compatible with non-zero winding num- 
ber. We see that skyrmion solutions can be interpreted 
in the above manner along with the additional condi- 
tion that u and v are mutually othogonal and satisfy 




FIG. 3: For an arbitrary smooth magnetization configura- 
tion, contour lines of n z and the phase of h x + ih y provide a 
natural coordinate system. For skyrmion solutions of the or- 
dinary ferromanget, this gives an orthogonal coordinate sys- 
tem with contour lines intersecting at right angles. Using 
the same ansatz for spinor condensate ferromagnets where 
contour lines intersect at right angles allows us to solve the 
non-linear and non-local equations of motion. The magneti- 
zation and contour lines are shown for the single skyrmion 
(left) for the ordinary ferromagnet and neutral configuration 
(right) for the spinor condensate ferromagnet. Hue indicates 
orientation of h x , h y components of the magnetization and 
brightness gives the h z component with white (black) indi- 
cating n z = +1 (n z = —1). 

Laplace's equation. Physically, these additional condi- 
tions on u and v can be understood as a consequence of 
minimizing the spin stiffness. The relationship between 
contour lines, winding number, and the exact solutions 
we discussed in this section is shown in Fig. |3j 

With this viewpoint, we can now generalize the 
skyrmion solutions for ordinary ferromagnets and also 
find new solutions for spinor condensate ferromagnets. 
For f(x + iy) — u(x, y) + iv(x, y) holomorphic, we take 

a = a(u), ft = kv (14) 

for the parameterization of Eq. [2] Compared to the 
skyrmion solutions of Eq. [3] with tan(a/2) = exp[u] for 
ordinary ferromagnets, we allow for general dependence 
a(u) for spinor ferromagnets. Whereas for ordinary fer- 
romagnets we only need to solve V 2 fi = 0, for spinor 
condensate ferromagnets we need to solve Eqs. [12l [13] 
with dt — 0. 

In addition, we include a constant of proportionality 
P = kv instead of P = v. Recall for ordinary ferromag- 
nets, P = +v and f(x + iy) holomorphic give skyrmion 
solutions with positive skyrmion density q while P = —v 
and f(x — iy) antiholomorphic give anti-skyrmion solu- 
tions with q. We can treat both types of solutions with 
just f(x + iy) holomorphic by allowing p — kv with k 
positive or negative. 

There are two cases to consider for f(z). The first is 
when f(z) is a polynomial in z with no singularities. This 
will turn out to describe skyrmion and anti-skyrmion 
stripe and domain wall configurations for both ordinary 
and spinor condensate ferromagnets. The second case is 
when f{z) has singularities. Since n should be single- 
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valued, j3 and thus kv can only have constant 2irN dis- 
continuities with N integer. This implies f(z) can only 
have logarithmic singularities. These solutions will turn 
out to simply be the localized skyrmion configurations for 
ordinary ferromagnets and neutral collections of localized 
skyrmions and anti-skyrmions for spinor condensate fer- 
romagnets. 

Next we consider Eq. [T2]for the superfluid velocity v. 
For the parameterization in Eqs. [TJ] we see from Eq. 3] 
that the skyrmion density q only depends on u. We thus 
take to only depend on u which reduces the equation 
— V 2 $ = q to — $"(it) = q(u). From here on primes 
denote derivatives with respect to u. By solving for <£>(ii) 
we can then obtain v by differentiating. Explicitly, we 
obtain 



nian in Eq. [TTJ is given by 



9 = -*cos(a)'|0,/| a , 



iSk[C + cos(a)]d z f (15) 



where w z = -v x — iv y , d z f = d x f — id y f and C is a con- 
stant of integration physically describing a u independent 
constant contribution to the superfluid velocity. 

We now proceed to analyze Eq. [T5] for spinor conden- 
sate ferromagnets. By substituting the results of Eq. [TBI 
above and the parameterization in Eqs. [21 1141 we find the 
z component of Eq. [12] is automatically satisfied. In ad- 
dition, the x and y components are proportional to each 
other and reduce to a second ordinary differential equa- 
tion for a(u). For completeness, we can use the same 
approach to analyze Eq. [TJ for ordinary ferromagnets us- 
ing the same parameterization in Eqs. [51 [3] 

For ordinary and spinor condensate ferromagnets, the 
equations of motion in Eq. [TJ and Eq. [T^] reduce to 

2a" = +k 2 sin(2a) 

2a" = -ASCk 2 sin(a) - (2S - l)k 2 sin(2a) (16) 

respectively. Notice the equations of motion for ordinary 
ferromagnets are formally given by the S = limit for 
spinor condensate ferromagnets. Recall the spin stiffness 
term scales linearly with S whereas the superfluid kinetic 
energy scales quadratically with S. For spinor condensate 
ferromagnets in the limit S — > 0, the superfluid kinetic 
energy is negligible compared to the spin stiffness and 
the ordinary ferromagnet is recovered. From here on, we 
consider the more general equation of motion for spinor 
condensate ferromagnets. 

Interpreting u as time, this equation is that of a clas- 
sical particle with coordinate a and momentum a' . The 
total energy E = K + U is a constant of motion where 
K = (a') 2 /2 is the kinetic energy while the periodic po- 
tential U and equation of motion are 

29—1 

U(a) = -2SCk 2 cos(a) — fc 2 cos(2a) 



x = y/E - U(a) 



(17) 



H = X I dudv^ 



2S+1 + 4C 2 S 



k 2 -E 



(18) 



with S = for ordinary ferromagnets. For this type of 
solution, the total energy calculated from the Hamilto- 



and similarly, S = for the term in brackets for ordinary 
ferromagnets. 

The integer N comes from changing variables (x, y) 
to (it, v) and taking into account each (u, v) may occur 
for multiple (x,y). Mathematically, it is given by the 
degree of / viewed as a map from the complex plane to 
itself C — > C. For localized skyrmion solutions for the 
ordinary ferromagnet, it physically corresponds to the 
skyrmion number. As an example, f(z) = \ogz for the 
solutions shown in Fig. [TJ(see Eq. [2lJfor a generalization) 
and each value of (u, v) occurs exactly once for (x, y) 
ranging over the plane and thus N = 1. For f(z) — 
log[(0(* - A, i)d{z - A*,i))/(0(* + A, i)tf(z + X*,i))} with 
$(z, t) the elliptic theta function and A = (l+i)/2 for the 
solutions shown in Fig. [5] (see Eq. [52]for a generalization) 
and each value of (u, v) occurs exactly twice for (x,y) 
ranging over the unit cell and thus N — 2. 

Here we comment on the significance of the parame- 
ters C and E. Different values of C and E correspond to 
solutions with different boundary conditions. Physically, 
C controls a constant background contribution to the su- 
perfluid velocity. The parameter E controls the relative 
scaling of the two components of the coordinate system. 
For example, for doubly periodic stripe solutions which 
we will discuss later, E controls the aspect ratio of the 
unit cell. 

Formally, E and C are constants of integration for the 
equations of motion. Since Eq. [T5] is a second order 
differential equation, we have to specify both a(u) and 
a'(u), with the latter given indirectly by the constant of 
motion E. In addition, C enters through integration of 
Eq. [15] relating the skyrmion density to the superfluid 
velocity. 

These considerations mean that the total energy in Eq. 
[T51 cannot be directly compared for different C and E. In 
particular, one should not consider minimizing the total 
energy % with respect to E and C. Specific values of E, 
and C will be selected by terms beyond the non-linear 
sigma model considered in this paper. 



A. Localized skyrmions and anti-skyrmions 

We begin by considering the case of localized 
skyrmions for ordinary ferromagnets and neutral collec- 
tions of localized skyrmions and anti-skyrmions for spinor 
condensate ferromagnets. This corresponds to f(z) hav- 
ing logarithmic singularities. The singularities should be 
of integer magnitude and k should also be an integer. 

Requiring a well-behaved, finite energy solution gives 
rise to several constraints. Consider Eq. [15] for the 
skyrmion density q and supefluid velocity v. Since 
d z f diverges at the logarithmic singularities, we require 
cos(a)' and C + cos(a) to vanish. In addition, a finite 
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region in (x, y) near logarithmic singularities is mapped 
to an infinite region in (u, v). The constant term in Eq. 
[TBI for the total energy is then integrated over an infinite 
interval. Thus, we also require it to vanish. 

These constraints uniquely specify E, C and the 
asymptotic value a-oo- For spinor condensate ferromag- 
nets, E = fc 2 (l + 65)/4, C = ±1, cos(a_oo) = =Fl with 
S = for ordinary ferromagnets. For C — ±1, we find 
for spinor condensate ferromagnets the solution of Eq. 
[TBI and the total energy of Eq. [TS] given by 



a(u) = cos _1 (=Ff) - 2cot~ 1 (v / 2S'sinh(fcu)) 
H = AirNSk 2 



1 



25tan- 1 (V25- 1) 



V25- 1 



(19) 



which describes either a — > 2ir or — 7r — > +7r kink solu- 
tion. For C — ±1 we find for ordinary ferromagnets 



a(u) = cos _1 (Tl) + 2tan- 1 (e fctl ) 
"H = 2nNSk 2 



(20) 







which describes in contrast either a — > tt or — 7r 
kink solution. 

Consider the classical mechanics problem describing 
the evolution of a(u) in Eq. [T7l For these solutions, 
E lies at the maximum giving rise to kink solutions. 
For ordinary ferromagnets, the kinks connect — > tt or 
— 7r — > and carry net positive or negative skyrmion 
charge, respectively. For spinor condensate ferromagnets, 
the kinks connect — tt — > +tt or — > 2tt and carry net 
neutral skyrmion charge. The neutral configurations con- 
sist of regions of oppositely charged skyrmion and anti- 
skyrmions. These regions are separated by lines where 
the skyrmion density q vanishes, the magnetization h is 
along z, and the superfluid velocity v is large. 

For f[z) having a finite number of logarithmic singu- 
larities, we can write 



f(z) = log 



Un=l( Z - a n) 



n 



N,, 



b m ) 



(21) 



with the degree N of the function f(z) given by N = 
m&x(N a , Nb). Here, a n (b n ) give the locations of +z 
(— z) merons each carrying net skyrmion charge 2tt for 
the ordinary ferromagnet. In contrast, a n (b n ) give the 
locations of skyrmions (anti-skyrmions) each carrying net 
skyrmion charge +4n (— 4tt) for the spinor condensate 
ferromagnet. We show the corresponding plots of n, q, 
and v in Fig. [JJfor f(z) — log(z). The single skyrmion 
solution for the ordinary ferromagnet with S — is shown 
on the left and a neutral configuration of one skyrmion 
and one anti-skyrmion for the spinor condensate ferro- 
magnet with S = 1 is shown on the right. 

For a periodic lattice of logarithmic singularities, 



/O) = log 



1X^1^(2 - On,r) 



/o 



n^=i^-6m,r) 



(22) 



5>V 



t 4 i v t «c 

• 4 4 % 



FIG. 4: Neutral stripe configuration satisfying non-trivial 
corner boundary conditions. Red (blue) background indicates 
positive (negative) skyrmion density q, black 2D arrows the 
superfluid velocity v, and shaded 3D arrows the magnetiza- 
tion h. 



where /o is a constant, and 1, r give the basis vectors 
generating the latttice in complex form, and t?(z, r) is 
the elliptic theta function. The elliptic theta function 
r) is essentially uniquely specified by the quasiperi- 
odic condition 

•d(z + n + niT, t) — exp(— irim 2 T — 2irimz)$(z, r) (23) 

and holomorphicity. Just as Eq. [21] is built up from the 
linear polynomials (z — zq) which are holomorphic and 
vanish at one point in the complex plane, Eq. [22] is built 
up from i?(z, r) which are holomorphic and vanish at one 
point in the unit cell. For a discussion of theta functions 
in the quantum Hall effect, see Ref. [43]. 

In the lattice case, N a = AT 6 and J2 a n — Yl b n hi 
order to have f(z) periodic. This restriction comes from 
requiring f(z + n + mr) — f{z) and using Eq. [521 The 
degree N of the function f(z) per unit cell is given by 
N = N a = Nb- Again, a n , b n give the locations of merons 
(skyrmions or anti-skyrmions) for the ordinary (spinor 
condensate) ferromagnet. Fig. [5] shows plots for f(z) 
having a lattice of logarithmic singularities with N — 
N a = N b = 2, f = 1, r = i, oi = -oa = (1 + i)/2, h = 
—62 = (1 — i)/2. Again the ordinary (spinor condensate) 
ferromagnet is on the left (right). 



B. Stripe configurations 

Now we turn to case of stripe configurations described 
by f(z) polynomial in z. The behavior of a(u) solutions 
controlled by the potential in Eq. [T7l changes as E crosses 
critical points dU /da = of the potential. 

For completeness, we first briefly consider f(z) given 
by higher order polynomials. The corresponding stripe 
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solutions are not doubly periodic, but satisfy non-trivial 
boundary conditions. For example, we show the mag- 
netization n, skyrmion density q, and superfluid veloc- 
ity v for f(z) = iz 2 in Fig. |U This solution satis- 
fies corner boundary conditions with zero normal compo- 
nent to both the superfluid velocity v and spin current 
J*, = nV p - e«*W'V ' li h k /2. 

From here on, we focus on f(z) = iz with the corre- 
sponding solutions doubly periodic and describing stripe 
configurations. The different types of behavior for a(u) 
are illustrated schematically along with the resulting con- 
figurations for n in Fig. [SJ E above the global maximum 
corresponds to a(u) monotonic in u which from here on 
we denote as M. This solution describes a periodic stripe 
solution with n z varying over the entire range ±1. E at 
a local maximum corresponds to a kink solution for a(u) 
connecting ot\ to ct2 denoted as K%*. This solution de- 
scribes a single domain wall configuration in n z and is 
the analog of the localized solution described earlier. E 
below a local maximum corresponds to a(u) oscillating 
near a fixed value ao denoted as O ao . Finally, E below 
the global minimum is forbidden, denoted as F. 

Notice that for 5 = 0, 5 = 1/2, 5 > 1/2, the cos(2a) 
term in the potential U(a) of Eq. [17] is negative, zero, 
and positive. For the montonic solutions M, this does not 
affect the qualitative behavior of the resulting periodic 
stripe configurations. For kink solutions K connecting 
— > 2tt or —it — > +7r (0 — > it or — tt — » 0), the result- 
ing single domain wall carrys zero net skyrmion charge 
(positive or negative skyrmion charge) for spinor con- 
densate (ordinary) ferromagnets. Oscillatory solutions 
O also have different behavior with oscillations centered 
about h z ss ±1 (h z w 0) for spinor condensate (ordinary) 
ferromagnets. 

For ordinary ferromagnet with 5 = 0, we parametrize 
E = fc 2 (l + 2<5)/4 and find the solution of Eq. [16] and the 
total energy of Eq. fTS] given by 

a(u) = tt/2 ± am (kuVTTS, (1 + Sy 1 ^ 



n 



2 



+ 0(1 + S) 



(24) 



where T-L is the total energy density given by the aver- 
aging % over the unit cell. Also, am(i, m) is the Jacobi 
amplitude function and we define 



9{m) = m 



Re^m" 1 )] 
Re^m" 1 )] 



(25) 



with K (E) the complete elliptic integral of the first (sec- 
ond) kind. For 5 < — 1 the solution is forbidden F. For 
— 1 < S < there are two oscillatory solutions at the 
±7r/2 minima 0± n /2- For 5 = there are two kink solu- 
tions Kg* and K ° w . For 5 > the solution is monotonic 
M. We show the classification of solutions along with the 
total energy density for ordinary ferromagnets in the bot- 
tom left of Fig. Notice the solutions and total energy 
density do not depend on C. This is because C enters 



through the superfluid velocity which is absent from the 
Lagrangian for ordinary ferromagnets. 

For spinor condensate ferromagnets with S = 1/2, we 
parameterize the constant of motion E = |C|fc 2 (l + 25) 
and find the solution of Eq. [16] and the total energy 
density from Eq. 1181 given by 



(u) = cob _1 (C/|C|) + 2am (kuy/\C\(l + S), (1 + 5)- 1 
(\C\-1) 2 5\C\ 



H 



\C\9(l + 6) 



(26) 



depends only on \C\. For 5 < — 1 the solution is forbidden 
F. For — 1 < 5 < there is one oscillatory solution O w 
(O ) for C < (C > 0). For 5 = there is one kink 
solution Kfi 71 (KZZ) for C < (C > 0). For 5 > the 
solution is monotonic M. We show the classification of 
solutions along with the total energy density for S = 1/2 
spinor condensate ferromagnets in the bottom right of 
Fig. [1 

For spinor condensate ferromagnets with S > 1/2 we 
parameterize the constant of motion E = (2S — l)k 2 (l + 
25) /A and the constant C = (25 - 1)7/25. With 7 = r 
and 5 = — 1 + 2ar where a = ±1 we find one solution for 
Eq. [16] while the total energy density from Eq. [18] given 
by 



a(u) = Arg 



t — a + 2\J t(t — a)s(u) + ts(uY 
I + <jt — <7ts(u) 2 

[h Q + (25 - l)jo] 



where we define the auxiliary function 



s(u) = sm(kuy/(l - ar)(2S - 1)) 
in the solution for a(u) and the functions 
((25- 1)t-ct5) 2 



(27) 



(28) 



h 



Jo 



45 

ot\Jt — a < (tt < 1 
otherwise 



(29) 



for the total energy density. With 7 = r sinh(r) and 
5 = —1 + 2r sinh(r) we find two solutions for Eq. 1161 and 
the total energy density from Eq. 1 181 given by 



a±(u) =Arg 



cosh(r/4)e 



=F sinh(r/4) 



sinh(r/4)e ±ra (") =F cosh(r/4) 
[h + (25 — l)j T (25 — l)fc] 



where we define the auxiliary function 

1 + r — r cosh(i) 



w(u) =am 



kus/2r{2S- 1), 



2r 



(30) 



(31) 
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FIG. 5: Stripe configurations for different boundary conditions in S > 1/2 spinor condensate (bottom row) and S = ordinary 
ferromagnets (top row). Red (blue) background indicates positive (negative) skyrmion density q, black 2D arrows the superfluid 
velocity v, and shaded 3D arrows the magnetization n. Notice h x , n y (n z ) wind horizontally (oscillate vertically). The first 
column shows schematic plots of the classical periodic potential U(a) controlling the evolution of the angular variable a for 
h z = cos(q). Labels for the corresonding type of solution are below the dashed lines indicating the corresponding constant of 
motion E. 

For E above the maximum of U, a(u) is montonic in the coordinate u giving rise to periodic stripe configurations labeled M 
with n z covering the entire range ±1. For E at a potential maximum, kink solutions connecting the maxima ati, ct2 give rise 
to single domain wall configurations in h z labeled . For a given value of C, kink solutions occur at one specific value of 
E. For E below the maximum, a oscillates about a minimum located at ao giving rise to periodic stripe configuration with h z 
oscillating about cos(ao). 

Notice that monotonic solutions M are qualitatively similar for both S > 1/2 spinor condensate and S — ordinary ferromag- 
nets. In contrast, spinor condensate (ordinary) ferromagnets have a single (two distinct) 2ir (ir) kink solutions K centered about 
a = (a — ±7r/2). In addition, there is a single (two distinct) oscillatory solutions O also centered about a = (a = ±7r/2) 
for spinor condensate (ordinary) ferromagnets. 



in the solutions for a(u) and the functions 

_ 1 - (2S- 1)6 (2.S- l) 2 /3 2 
1 ~ 2 + AS 

. _ Re[2V?^(a;)] 

Re[2VF[cosh(V2) 2 n(-sinh(i/2) 2 ,w) - K(u) 

, J ^rWSTV O„O Phase 
I otherwise 



1 



r cosh(r) 



2r 



(32) 



for the total energy density where IT(m, n) is the complete 
elliptic integral of the third kind. We show the classifica- 
tion of solutions along with the total energy density for 
S > 1/2 spinor condensate ferromagnets in the top row 
of Fig. [B] The boundaries between solutions of different 
types are given by 5 = — 1 + 2/3, S = -1 - 2/3, S = /3 2 . 



For increasing 7, notice kink solutions evolve from just 
one Kq 7 ' through a region with two K^~ a and K*%> 
to just one K±l for 7 < -1, -1 < 7 < +1, +1 < 7, 
respectively. The kink solutions separate the monotonic 
solutions M from the oscillatory solutions. For increas- 
ing 7, the oscillatory solutions also evolve from just one 
to a region with two O v and Oq, to just one Oo for 
S < -1 - 2 7 , -1 + 2 7 > 6 > -1 - 27, S < -1 + 2 7 , 
respectively. 



VI. DISCUSSION 

Having presented a unified description of both local- 
ized and extended stripe solutions in ordinary and spinor 
condensate ferromagnets, we now turn to how these solu- 
tions offer insight into different physical phenomena. We 
first consider quantum Hall systems. As discussed in Sec. 
IIII1 configurations for quantum Hall ferromagnets away 
from quantum Hall plateaus carry net skymrion charge 
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FIG. 6: Classification of solutions (2D plots) and total energy density (3D plots) for stripe configurations given different 
boundary conditions in S > 1/2 (top row) and S — 1/2 (bottom left) spinor condensate as well as S = (bottom right) 
ordinary ferromagnets. Fig. [S] illustrates the corresponding configurations. 7 controls a constant contribution to the superfluid 
velocity. 5 controls the energy of an associated classical mechanics problem giving the evolution of h z = cos (a), the z component 
of the magnetization. Monotonic solutions M have h z covering the entire range ±1. Kink solutions K"^ describe single domain 
wall configurations with n z connecting cos(qi) to cos(a2). Oscillatory solutions O aQ have n z oscillating about cos(«o). Notice 
kink solutions K always separate monotonic M from oscillatory O solutions. For S > 1/2 spinor condensate ferromagnets, 
notice the two distinct oscillatory and kink solutions for each 7, 5 in the region near the origin which are absent for S = 1/2. 
For S = ordinary ferromagnets, there is no dependence on C. 



[HIM USUI, S3- Thus > solutions for the ordinary fcrro- 
magnet describing collections of localized skyrmions car- 
rying net charge as shown in the left of Figs. Q] and [5] 
have been used extensively in this regime. 

However, we showed in Sec. |V]that these solutions of 
localized topological objects can be derived in a unified 
framework along with extended stripe solutions. There 
have been a number of studies on the possibility of quan- 
tum Hall states with stripe order. At high Landau levels 
and with frozen spin degrees of freedom, Coloumb inter- 
action may directly favor charge density waves as pre- 
dicted theoretically 0, HI| and verified experimentally 
[iti H3] . Such states are not directly comparable to the 
stripe solutions we describe which have fixed total den- 
sity and stripe order in the relative density. However, 
stripe order has also been proposed [48l450|| and experi- 
mental evidence observed [51[ in the context of quantum 
Hall bilayers. Here, even though the total density be- 
tween layers is fixed, both interlayer coherence and rela- 
tive density imbalance can develop. The isospin degree of 
freedom that arises can be used to define an appropriate 
magnetization vector h. Here, the phase of the interlayer 
coherence gives the orientation of h x , n y , while the rel- 
ative density imbalance gives rl z . States with skyrmion 



stripe order and winding h x , n y have been proposed that 
are direct analogs of the configurations shown in the top 
row of Fig. [5j 

For spinor condensate ferromagnets, experiments at 
Berkeley suggest the possibility of a condensate with 
crystalline magnetic order [l5|, [l6[ . This crystalline or- 
der arises from an effective dipolar interactions modified 
by rapid Larmor precession and reduced dimensional- 
ity. It can drive dynamical instabilities of the uniform 
state which occur in a characteristic pattern [lj| fl8j . 
Modes controlling the component n parallel (perpendic- 
ular) to the magnetic field in spin space are unstable 
along wavevectors perpendicular (parallel) to the mag- 
netic field in real space. Instabilities of this type can give 
rise to the spin textures shown in the stripe solutions of 
the bottom row in Fig. [SJ Here, h z is modulated along 
the y direction while h x , h y wind along the x direction. 

In the companion paper |21| . we have performed a 
systematic numerical study of minimal energy configu- 
rations for spinor condensate ferromagnets with dipolar 
interactions. This is made possible by the use of sym- 
metry operations combining real space and spin space 
operations to distinguish different symmetry classes of 
solutions. For applied magnetic field in the plane B = x 
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FIG. 7: Numerically optimized configuration for two- 
dimensional spinor condensate ferromagnets with an effective 
dipolar interaction modified by rapid Larmor precession. The 
magnetic field B = x inducing Larmor precession lies along 
the horizontal axis in the plane. Lattice constants are ay = 90 
fim and a± — 42 /im. Red (blue) background indicates pos- 
itive (negative) skyrmion density q, black 2D arrows the su- 
perfluid velocity v, and shaded 3D arrows the magnetization 



corresponding to current experiments, we show the low- 
est energy configuration in Fig. [7] 

Notice n z is modulated between ±1 just as in the 
monotonic M solutions shown in Fig. [5] In addition, 
the n x , fiy components wind along the horizontal axis. 
However, notice the winding in the h x , h y components is 
not uniform as in the solutions we find in this paper. In 
addition, the winding changes from clockwise to counter- 
clockwise halfway along the horizontal axis. For the so- 
lutions we find in this paper, the skyrmion density forms 
stripes of opposite charge parallel to the horizontal axis. 
For the lowest energy configuration in Fig. the non- 
uniform winding leads to concentration of the skyrmion 
density in smaller regions and modulation in the sign of 
the skyrmion density along the x axis. We find minimal 
energy configurations in other symmetry classes are gen- 
erally of this type with n z oscillating between ±1 along y 
and h x , n y winding along x. However, the detailed form 



of the winding along x varies for different classes. 

Thus we see that the exact solutions for spinor conden- 
sate ferromagnets without dipolar interactions provides 
a more transparent physical picture for the numerical so- 
lutions with dipolar interactions. This can be seen as fol- 
lows. The solutions we find in this paper describe overall 
neutral collections of skyrmion and anti-skyrmion topo- 
logical objects. The neutrality constraint comes from the 
long-ranged divergence of the skyrmion interaction and 
remains even when considering additional spin interac- 
tions such as the dipolar interaction. Morever, skyrmions 
and anti-skyrmions themselves have a non-trivial spin 
texture which is evident in Fig. [5]showing the solutions of 
this paper. When dipolar interactions are included, such 
spin textures can take advantage of the gain in dipolar in- 
teraction energy without chaning their qualitative struc- 
ture. However, quantitative details for minimal energy 
configurations such as the one shown in Fig. [7] require 
detailed analysis of the competition between dipolar in- 
teractions, skyrmion interactions, and spin stiffness. 

In conclusion, we have presented the low-energy effec- 
tive theory of spinor condensate ferromagnets. This ef- 
fective theory describes the superfluid velocity and mag- 
netization degrees of freedom and can be written as a 
non-linear sigma model with long-ranged interactions be- 
tween skyrmions, the topological objects of the theory. 
Quantum Hall ferromagnets share a similar effective the- 
ory with long-ranged skyrmion interactions. For the case 
of spinor condensate ferromagnets, we find exact solu- 
tions for the non-linear equations of motion describing 
neutral configurations of skyrmions and anti- skyrmions 
carrying zero net skyrmion charge. These solutions de- 
scribe within a unified framework both collections of lo- 
calized topological objects as well as extended stripe con- 
figurations. In particular, they can be used to understand 
aspects of non-trivial spin textures in both quantum Hall 
ferromagnets as well as spinor condensate ferromagnets 
with dipolar interactions. 
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